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Abstract 

We consider the random matrix model with external source, in case where the poten- 
tial V(x) is an even polynomial and the external source has two eigenvalues ±a of equal 
multiplicity. We show that the limiting mean eigenvalue distribution of this model can be 
characterized as the first component of a pair of measures (fii, ^2) that solve a constrained 
vector equilibrium problem. The proof is based on the steepest descent analysis of the 
associated Riemann-Hilbert problem for multiple orthogonal polynomials. 

We illustrate our results in detail for the case of a quartic double well potential V(x) = 
jX 4 — \x 2 . We are able to determine the precise location of the phase transitions in the 
ta-plane, where either the constraint becomes active, or the two intervals in the support 
come together (or both). 



1 Introduction 

The random matrix model with external source is the probability measure 

— exp (— n Tr(V(M) — AM)) dM (1.1) 

defined on the space ofnxn Hermitian matrices M. Here A is a given Hermitian matrix (the 
external source), V : R — >■ R is a function with sufficient increase at infinity (the potential), and 
Z n is a normalization constant. 

The model (|1.1[) was first studied by Brezin and Hikami jTOJ [TT] and P. Zinn- Justin [55J IHH] 
who showed that the eigenvalue correlations are determinantal. In [7] it was observed that the 
correlation kernel can be expressed in terms of multiple orthogonal polynomials. Due to the 
Riemann-Hilbert problem for multiple orthogonal polynomials |37) this opened up a new way for 
asymptotic analysis. For the quadratic case 

V(x) = \x 2 (1.2) 

with external source 

A = diag(a, . . . , a, —a, . . . , —a) (1-3) 

n/2 times n/2 times 
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with two eigenvalues ±a of equal multiplicity (thus n is even), this was done in great detail in 
the three papers [1IE1IH]- 

The quadratic case is of special interest because it has an equivalent formulation in terms of 
non-intersecting Brownian motions that start at one value and end at certain prescribed values 
[2] which is a variation on Dyson's Brownian motion |22j . The quadratic model with external 
source (|1.3[) exhibits a phase transition, since for small a > the eigenvalues accumulate on one 
interval while for larger a the eigenvalues accumulate on two disjoint intervals. At the critical 
value of a the local eigenvalue correlations are given in terms of Pearcey integrals [TJ ESI EH] • 

In this paper we study the external source model (jl.ip with a more general potential V. We 
assume that V is an even polynomial 

d 

V{x) = ^v jX 2 \ v d >0 (1.4) 
j=i 

of degree 2d. 

For a = the external source model (11.11) reduces to the usual unitary matrix model 

^-exp(-nTrV(M)), (1.5) 
An 

which is one of the most studied models in random matrix theory in both mathematics and 
physics, see e.g. [5J [T31 US] for rigorous study using the Riemann-Hilbert approach. A 
basic fact is that for n — > oo, the limiting mean eigenvalue distribution of the matrix M in (|1.5[) 
minimizes the energy functional 

E(n) =U log -±-dfjL(x)dn{y) + J V(x)dn(x) (1.6) 

over all probability measures fi on R. 

It is an open problem to find an analogue for the equilibrium problem (I1.6P in the general 
context of the random matrix model with external source. This paper contains a first result 
in this direction. We consider the external source model (jl.ll) in case where the potential V is 
an even polynomial (|1.4p . The external source A is again given by (jl.3D with two eigenvalues 
±a of equal multiplicity. We show that under these assumptions, the limiting mean eigenvalue 
distribution of the matrix M in (jl.ll) exists, and that it arises as the first component of a pair 
of measures (/Xi,^) solving a certain vector equilibrium problem, see Section [2] 

We will illustrate our results in detail for a particular case of a non-convex potential, namely 
the quartic double well potential 

V{x) = - A x i - t -x 2 , t>0. (1.7) 

For the quartic model (|1 .5(1 . fj 1 . T[) (without external source) it is known that the eigenvalues 
accumulate on either one or two intervals. The local eigenvalue correlations for the critical value 
of t = t cr = 2 are given in terms of ^-functions associated with the Hastings-McLeod solution 
of the Painleve II equation [51 [T^] . 

So in the quartic model with external source there exist at least two mechanisms by which 
a transition from one to two intervals can occur: namely a Pearcey transition and a Painleve II 
transition. It will be one of the outcomes of the present paper that we can determine precisely 
the location of the phase transitions in the ia-planc. 
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2 Statement of results 



2.1 Equilibrium problem 

The main ingredient in our analysis is a new vector equilibrium problem associated with the 
random matrix model Ql.lj) with external source. We emphasize that it only applies in the 
setting we are considering, namely an even polynomial potential V as in fj 1 -4[) and an external 
source (|1 .3|) with two eigenvalues of equal multiplicity. This setting gives a symmetry with 
respect to the origin, which we use in an essential way. 

The equilibrium problem is as follows. We minimize the energy functional 



E(fJ,x,H2) = JJ log l _ din[x)dn\{y) + JJ log , ^ , dm 



(x)dfj, 2 (y) 



log- — - — -d^i{x)d^ 2 {y) + [ (V(x) - a\x\)dfii(x) (2.1) 
\x-y\ J 

with respect to all pairs of measures (/ii, ^2) satisfying 

• fi\ and fi2 have finite logarithmic energy, 

• fi\ is a measure on R with total mass 1, 

• ^2 is a measure on iM. with total mass 1/2 that satisfies the constraint 

M2 < a (2.2) 
where a is the measure on iM with constant density 

dcr a _ . . 

rn = - «e«. 2.3 



Standard references on potential theory in the complex plane are [MJ [321 [35] . 

The equilibrium problem (|2.1[) has both an external field V(x) — a|a;| acting on fix, and an 
upper constraint a acting on The interaction between fi± and /12 is of Nikishin type [26] . 
This type of vector equilibrium problem also appeared recently in a model of non-intersecting 
squared Bessel paths [5S] and in the two-matrix model with quartic potential [2"0] . 

Our first result concerns the structure of the minimizer of the equilibrium problem. 

Theorem 2.1. There is a unique minimizer (//i, ^2) which satisfies 
(a) The support of /ii is bounded and consists of a finite union of intervals 

N 

sapp(jii)=\J[a j ,b j ]. (2.4) 

3=1 

The measure /zi is absolutely continuous with density 
dpi (x) 



dx 



h(x) 



N N 



Y[(b j -x)(x-a j ), x€\J[aj,bj], (2.5) 
\ 3=1 3=1 



where h is a nonnegative function on supp(/xi) = [JjLi i a j'bj] that is real analytic, except 
possibly at zero. 
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(b) The support of \ii is the full imaginary axis and there exists c > such that 

supp(er — jjbi) = {— ioo, —ic] U [ic, ioo). (2-6) 
We have that c = if and only if 

and in that case, /i 2 has the density 

d^ 2 {z) _ 1 
\dz\ ~ 2tt 

If (12.71) is not satisfied then c > is determined by the condition 

d/xi (s) 



d/ii(s), zeii (2.8) 



and in that case 



and 



y/s 2 + c 

d/it 2 (z) a 
\dz\ 7r' 



2a, (2.9) 



if z G [~ic,ic], (2.10) 



d M2 (z)_a 1 f \zW\z\*-<* d ^ (s) 



\dz\ it 2tt J (|z| 2 + s 2 )Vs 2 +c 2 



2nJ (| z |2 + s 2)^ 

if z G (—ioo, — ic] U [ic, ioo). 

(c) Both fj,i and /i 2 are symmetric with respect to the origin. 

Theorem [O will be proved in Sections O and EH Note that (|2"~j ]) -(|2lji |) are similar to 
statements proved in [20], while (|2.9[) and (|2.1ip - (|2.12|) have apparently not been stated before. 

2.2 Variational conditions 

The minimizcr (pi, /x 2 ) to the equilibrium problem in Scction ^.ll is characterized by the following 
Euler-Lagrange variational conditions. We write 



U»(x)= [ log y^—My) 
J F-J/l 



for the logarithmic potential of a measure /x. 

Proposition 2.2. The measures [i\ and /i 2 satisfy for some constant <€l, 

2U>* 1 (x) = U^ 2 {x)-V(x)+a\x\-e, x G supp(/ii), (2.13) 

2U tll (x) > U^{x) - V(x) +a\x\ - £, \ supp(^i), (2.14) 

2U tl2 (x) = U^{x), x G supp(cr-^ 2 ), (2.15) 

2U" 2 (x) < U^ix), xGiR\supp(CT-^ 2 ). (2.16) 

These relations follow directly from the variational conditions of the equilibrium problem. 
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2.3 Regular and singular cases 

We say that \i\ is regular if in (|2.5p we have h(x) > on supp(/ii) = ,bj], and if the 

variational inequality (|2.14p is strict for every x € K \ supp(/ii). Otherwise n\ is called singular 

maim. 

The measure fi2 has a density on iM. which is bounded by ~. We say that \i2 is singular if 
equality in this restriction is attained at z = and at no other point of iM. In all other cases (in 
particular if c > 0), the measure fi 2 is called regular. 

For the analysis in this paper, we will assume that both //i and /i2 are regular. The following 
lemma follows immediately from this assumption and is stated only for further reference. 

Lemma 2.3. The measures [i\ and a — fi2 satisfy the following square root behavior near their 
endpoints dj, bj and ±ic 

(a) If the measure is regular then it has a density of the form (|2.5p with h strictly positive 
on \jf =1 [aj,bj]. 

(b) If c > then the measure a — \i% has a density of the form 

a d(i 2 (z) 



where k is an analytic, strictly positive function on (— ioo, —ic] U [ic, ioo). 

Lemma 12.3( a) follows immediately from the definition of ji\ being regular. Lemma 12.3( b) 
follows from (|2.11l) . 



2.4 Limiting eigenvalue distribution 

Our main result deals with the global distribution of eigenvalues as n — > oo. 

Theorem 2.4. Let V be an even polynomial, and let A be a diagonal matrix with two eigenvalues 
±a of equal multiplicity. Let{)i\,\i2) be the solution of the equilibrium vroblem in Section \2.1l and 
assume that both [i\ and /j,2 are regular in the sense explained above. Then the mean eigenvalue 
distribution of a matrix M from the random matrix model 

— exp (—n Tr(V(M) - AM)) dM 

has the limit fi\ as n — > oo. 

We strongly expect that the conclusion of Theorem 12.41 remains valid in the case where /zi 
and/or \i2 is singular. 

Theorem 12.41 will be proved in Section 15771 



2.5 About the proof 

The proof of Theorem 12.41 is based on the Riemann-Hilbert problem for multiple orthogonal 
polynomials and its connection with the external source model (jl.lj) . 

The multiple orthogonal polynomials P n (x) in question are orthogonal with respect to the 
weights 

tui(a:) = e - n ^- ax \ w 2 {x) = e-< v ^ +ax \ (2.17) 
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More precisely, P n (x) is a monic polynomial of degree n that is characterized by the multiple 
orthogonality conditions (we assume n is even) 

P n {x)w j {x)x k dx = 0, j = l,2, fc = 0,...,n/2- 1. (2.18) 

The polynomial P n (%) is also the average characteristic polynomial 

P n («) = E [det(zl n - AT)] , (2.19) 

where the expectation E is taken with respect to the model see [7J. 

The Riemann-Hilbert problem (RH problem) for P n (x) is as follows. We look for an analytic 
matrix-valued function Y : C \ R — > C 3x3 with jump properties 

(1 wi(x) u> 2 (x) N 
1 | , x e E, (2.20) 

1 

and asymptotic condition 

y(z) = (J + 0(1/ z)) I z-"/ 2 I asz^oo. (2.21) 




Here and below we use the following standard notation. If T is an oriented contour in 
the complex plane, then the side that is on the left (right) when traversing T according to its 
orientation is called the +- side ( — side), and for any z € T we use Y+(z) (Y_(z)) to denote the 
limiting values of Y(z) along the +-sidc ( — side) of V. In (|2.20l) we have the contour r = K 
oriented from left to right. 

The RH problem fl2T2T)]) - (l2~2"T]) has a unique solution. The (1, l)-entry of Y(z) is the multiple 
orthogonal polynomial P n (z) characterized by (|2.18D . 

It is known that the eigenvalues of the random matrix model with external source (|1.3[) form 
a determinantal point process with correlation kernel [7J 

K n (x,y) = 2ni( l_ y) (0 My) My)) Y+^Y^x) . (2.22) 

Theorem 12.41 then comes down to the following statement about the limiting behavior of the 
kernels K n : 

lim -K n (x,x) = d ^ X \ (2.23) 
n-s-oo n dx 

We will establish (|2.23[) in Section [5771 thereby proving Theorem 12.41 

From the RH analysis it is possible to obtain universality results for the local eigenvalue 
correlations as well. In the regular cases, this leads to the usual sine kernel in the bulk and Airy 
kernel at the edge points of the spectrum. We will not discuss this any further and refer to the 
papers [3J [5J HI [TCJ [TTJ , among others, for a detailed analysis in a similar context. 

2.6 Organization of the paper 

The rest of the paper is organized as follows. In Section [3] we discuss the structure of the 
equilibrium measures and we prove Theorem 12. II In Section U we introduce the Riemann surface 



G 



built from the solution of the equilibrium problem. Section [5] contains the steepest descent 
analysis of the RH problem for Y(z), leading to the proof of Theorem l2.4l In Section |5] we make 
some general remarks on the expected phase transitions of our model, and in Section [7] we study 
this in detail for the case of a quartic potential. 



3 The equilibrium problem 

3.1 Existence of the minimizer 

In this section we prove the existence of the minimizer (hiiP-2) of the equilibrium problem in 
Section O To this end we follow [20j Section 4]. 



Proof. The energy functional (I2.1|) can be written as 

E(jix, 1*2) = I /(mi) + jAMi - 2/i 2 ) + / (V(x)-a\x\)d f x 1 (x) (3.1) 

where 

I{v) =11 log -dv{x)dv{y) 

J J \x-y\ 

denotes the logarithmic energy of a signed measure v. Occasionally we will also write 

I{vi,v 2 ) = ([ log- — - — -dv 1 (x)dv 2 {y) 
J J \x-y\ 

to denote the mixed energy of a pair of measures v\ and v 2 ■ 

Since I{y) > if v is a signed measure with J dv = 0, we find from (|3 . 1 [) that 

E(pi,H2)>~IQii) + J (V(x)-a\x\)dm(x) 



> min 

pi 



> — CO, 



where the last inequality follows from standard logarithmic potential theory with external fields, 
see e.g. |35j . Thus the energy functional is bounded from below. 

If we fix H2 on iR and we minimize (|2.ip with respect to fii only, then the problem for fi\ is 
to minimize 



/(Mi)+ / (V(x) - a\x\ - U^(x)) d^(x) 

The extra term U^ 2 (x) comes from the interaction between /ii and fi%. It is a term that attracts 
the Hi mass towards the origin. It can indeed be proved (as in [10]) that if the minimizer in 
external field V(x) — a\x\ is contained in [— X, X], then the minimizer in external field V(x) — 
a\x\ — C/ M2 (x) is also contained in [— X, X] (and so X is independent of (12)- 
If we fix [ii on [—X, X] then the problem for [12 is to minimize 

/(M2)- ju^{x)dH2{x) 

among all [12 < o~ with total mass 1/2. As in |20j . equality in the constraint is attained precisely 
on an interval of the form [~ic,ic] for certain c > 0. We will show further that the minimizer 
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\i2 satisfying this constraint is given explicitly by (|2.8I) - (|2.12I) . From these explicit formulas it 
follows immediately that for a measure [i\ on [— X, X], the corresponding minimizer /i2 satisfies 

du,2(z) K _ 

with a constant K that only depends on X . 

Let (/ii, n , /Z2,n)5£Li be a sequence of vectors of measures satisfying supp(//i, n ) C R, supp(/i2,n) C 
tR, J d/ii.n = 1, J d/i2,n == 1/2, and jLt2,n < 0", so that 

^(Ml,n,M2,n) < - +inf E(fJ, 1 ,fJ, 2 ). 

n 

As shown above, we may assume in addition that 

supp(/ii,„) C [~X,X] 
d^A z ) < K T( . m 

with X and K independent of n. 

Then it follows as in [50] that the sequences (Ati,n)^i and (M2,n)5?Li are tight. There is a 
convergent subsequence of (jUi^j M2,n 

and the limit is the vector of minimizing measures, see 

also [27] . 

Summarizing, we have now proved the existence of the solution (pi, /X2) to the equilibrium 
problem. The uniqueness of the solution follows in a standard way from the convexity of the 
energy functional, see e.g. (|3.ip and [5U]. □ 



3.2 Proof of Theorem EH] 

In this section we prove Theorem 12.11 

Proof. The proof of Theorem ^.lf a) follows as in |15j , while Part (c) is evident from the symmetry 
of the problem. 

It remains to prove (12.7[) - (|2.12[) in Theorem l2.1f b'). For k = 1, 2, define the Cauchy transforms 

Fk(z)= / dfik(s), z G C\supp(// fe ). (3.2) 

J z - s 

By differentiating the variational condition (]2 . 15[> we find that 

F 2+ (z) + F 2 -{z) - Fi{z) = 0, z G (-ioo,-ic) U (ic,ioo). (3.3) 

Here we assume that the imaginary axis is oriented from bottom to top, so that the +-side is 
on the left, and the — side is on the right, as usual. If c > then = ^tjt — — on (—ic.ic), 

' 1 \dz\ \dz\ 7T v ' 

from which it follows that 

F 2 +{z) - F 2 -{z) = -2a, z G (-ic,ic). (3.4) 

We can solve equations (|3.3|) — (|3.4|) for F2. We consider the two cases: c = and c > 0. 
Case 1, c = 0. In this case we have equation (|3. 31) on iR, 

F 2+ (z) + F 2 _(z)-F 1 (z) = 0, *G (-too, too). (3.5) 
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Define the Cauchy transforms of the restrictions of the measure [i\ to the positive and negative 
half-axes, 

+ = r°° Ws)^ zec\(R + n supper)) , 

Jo z-s 

Ff(z)= / zeC\(M_nsupp( Ml )). 

J — oo z s 

Then 

F 1 (z) = F+(z) + F 1 -(z), 
and, due to the uniqueness of the solution of the scalar Riemann-Hilbert problem (|3.5[) , we have 



F2(z) = 



F+(z), Rez<0, 
Ff(z), Rez>0. 



Therefore, the Sokhotski-Plemelj formula implies that for z G iK, 

' //( - >( : ' [ [F 2+ (z) F 2 _(z)] = -i- - Ff(z)] 



i r sjMi(>) 

izi 2 + S 2' 



i 

~~ ~2^ 

which is equivalent to (|2.8I) . The case c = is valid if and only if the density (12.81) is bounded 
by a/ ir. Since (|2.8|) assumes its maximum for |z| =0, this happens if and only if (|2.7[) holds. 
Case 2, c > 0. Consider the function 

F 2 (z) = ^jfi= i?(z) = z 2 + c 2 , 



where ^/ R(z) = y/ z 2 + c 2 is defined with a cut (— ioo,— ic] U [ic, ioo), and ^/ R(0) = c. From 
equation p. 31) we obtain that 

F 2+ (z)-F 2 _(z) = , ^ e (-ioo,-ic]U [ic,ioo), (3.6) 

V i? ( z )+ 

where + again denotes the limiting value from the left half plane. Observe that for y > c, 



\/R{iy) + = -Wy 2 - c 2 

In addition, from equation Q3.4[) we obtain that 



F 2+ (z) - F 2 -{z) = ze(-ic,ic). (3.7) 

By (|3.6[) and (|3.7[) the Sokhotski-Plemelj formula implies that 

F(z)^-^f r\ d Fi{s)ds i 1 r 2ads f 38 ) 
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m r 



Figure 1: The contour T. 



The first term in the right-hand side of (|3.8p can be written as 

1 ( r ic r°°\ Fi(a)da 1 f F 1 (s)ds 



2m 



— IOC J IC 



(3.9) 



where the contour T is depicted in Figure [T] From (|3.2I) and Fubini's theorem, 



Fi (s) ds 



ds 



1 

J r (z - s)y/WjB) 4mJ-ooJr (z - s)(s - Qy/R(s) 

i r°° i / r ds 



d M i(C) 



ds 



4:^^iJ_ 00 z-C\J^{z-s)^/R(s} J r (s - Qy/R(s) , 
By contour deformation and Cauchy's theorem we have that 

' z G C \ ((— ico, — ic\ U [ic, ioo)), 



d Ml (C). (3.10) 



47ri/ r (2_ fl ) A /fl(i) 2 v / i?(^)' 
hence (|3.10p reduces to 



Fx (s) ds 



47ri/ r (z-s)vARp) 



1 



1 



i /■<» oMO 



y/W)J-oc 2(«-C) J-oo 2(z-C)/R(C)' 

For the second term in the right-hand side of (|3.8I) we have 

1 /" lc 2a ds a sgn Re z 



2niJ_ ic ( z -s)^/R(s) y/R(z) 

By inserting (13.111) and f|3 . 1 2[) in (|3.8[) . we obtain that 

~ . . 1 f°° dfj,i(s) f°° d/zi(s) a sgn Re z 

F 2 {z) - 



(3.11) 



(3.12) 



V^i-oo 2(z-s) 7-oo 2(z-s)/R(Y) ' 
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tl6IlC6 

, , FUz) j-z f°° duds) 

Note that by taking z — > +00 in (|3.13ll . we find the relation (|2.9|) between c and a. Now the 
density of [12 is equal to 

d^2{z) 1 r , , , 
^ = -2^ [F2+(z) - F2 - (z)] 

_ a z^|z| 2 -c? f°° dm(s) 



, 2(z - s)Vs 2 + c 2 
a ^|z| 2 -c 2 f°° ( 1 1 \ d/il(«) 



7T 7T 7 \ Z - S Z + S J 2V S 2 + C 2 



a 1 p |z|^|z| 2 -c 2 c^i(s) 



7T 7T 



M 2 + « 2 V^Td^ 

which is equivalent to ([2XT]t . Then (|2J2|) follows from this and (|2J|) . □ 

3.3 Structure of the equilibrium measures in the regular case 

Theorem 12.11 implies that in the regular case, the structure of the equilibrium measures near 
the origin is described by one of the following three cases. This distinction will be important at 
several places of our RH steepest descent analysis. 

Case T. c — and N is even. 

Case II: c > and N is even. 

Case III: c > and N is odd. 

Here we recall the definitions of c and N in (|2.4[) - ([2.6[) . Thus Case I could be formulated 
cquivalcntly as supp(er — ^12) = and ^ supp(//i). Similar formulations can be given for 
Case II and Case III. Note that the situation where c = and N is odd is not in the above list 
since it would imply that u\ has singular behavior at the origin. 

For the quadratic potential V(x) = x 2 /2, it turns out that we are in Case I (with N = 2) for 
large values of a and in Case III (with N = 1) for small values of a. The Case II does not occur. 

In general, one expects Case I to happen for large values of a. Note that the upper constraint 
a is not active in this case and therefore it could be removed from the equilibrium problem. A 
consequence of this is that 2u2 is equal to the balayage of /ii onto the imaginary axis, see also 
HH5D . This means that 2U^-{z) = U^ 1 (z) for z e iR. 

Let Hi and fi^ denote the restrictions of /ii to the negative and positive real axis, respectively. 
By symmetry, we then have in Case I that U2 is the balayage of either /if or uf' onto the imaginary 
axis, and moreover 

C/ M2 (z) = U^(z), forRez>0, (3.14) 
U tl2 {z) = U^{z), forRez<0. (3.15) 

Note that the equality is valid not only on the imaginary axis, but also in a full half-plane. This 
follows from an easy application of the minimum and maximum principles for harmonic functions 
[33 Chapter 0]. 
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Then the following string of equations is easy to verify: 

= J U»?dni=I(tf,l*i), (3-16) 

and in a similar way 

/(/xi, /x 2 ) = / U^dtf + J U>»dn? = 2/(/if , /if ). (3.17) 

Since also 

/(/xx) = /(/if + /if) = /(/if) + /(/if) + 2/(/if , /if), (3.18) 
it then follows from (|3.16p - (|3.18p that the energy functional (|2.1I) can be rewritten in Case I as 

E(pL U /i 2 ) = J(Mi ) + /(/if) + /(Mi , ) 

+ / (y{x)+ax)dp%{x)+ I (V(x) - ax)dnf{x). 



Therefore the equilibrium problem is equivalent to the following equilibrium problem of Angelesco 
type for /if and /if: Minimize 

S(/if,Mf) = J J logj^-^dtf(x)dtf(y) + J J \ogj-L-^d(4(x)drf(y) 

log——dri( x )dfx?(y) 

\x-y\ 



+ J (V(x)+ax)dtf(x) + J (V(x) - ax)dfi? (x), (3.19) 

with respect to all pairs of measures (/if ,/xf) satisfying 

• /if is a measure on R_ with total mass 1/2 

• /if is a measure on R + with total mass 1/2. 

For the quadratic potential V(x) = \x 2 , this equilibrium problem was described in [5]. 

4 Riemann surface 

From the minimizer (/zi,//2) of the vector equilibrium problem we construct a three sheeted 
Riemann surface TZ, whose three sheets are given as follows. 

N 



TZ 1 =C\\J[a j ,b j ], (4.1) 
K 2 = C\ I [M+n {J[a jt bj] ] U{-ic,ic) ] , (4.2) 



j'=i 

JV 



ft 3 = c \ ( ( M_ n (J [o,-, &,] ] U (-tc, ic) ] . (4.3) 
i=i 
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Figure 2: Riemann surface for Case I (with N = 4). 




Here C = C U {00} denotes the Riemann sphere. 

The sheet TZi is connected with 72. 2 via the intervals [a 3 ,& 3 ] on the positive real line, IZi is 
connected with 72-3 via the intervals [a 3 ,b 3 ] on the negative real line, and (in Case II and Case III) 
IZ2 is connected to H3 via the interval [—ic, ic] on the imaginary axis. The connections are in 
the usual crosswise manner. The Riemann surface is compact and has genus 

fTV-2, in Case I, , . 

y \ N - 1, in Case II and Case III. K 1 

Here the Cases I, II and III were defined in Section |3~51 An illustration of the Riemann surface 
for each of these three cases is shown in Figures [2HU 

Recall the functions Fx and F2 in (|3.2p . These functions are used to define a meromorphic 
function on the Riemann surface, compare with [2Q[ Lemma 5.1]: 
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Figure 4: Riemann surface for Case III (with N = 3). 



Proposition 4.1. For k — 1,2, 3, /ei £fc(z) on the sheet TZk be defined by 



6(* 

£2(2 



= y'(z)-Fi(z), zeKi, 

= a + F^z) -F 2 (z), zeK 2 , Rez>0, 

= a + F 2 (z), zelZ 2 , Rez<0, 

= -a + F 2 (z), zeTZ 3 , Rez>0, 

= -a + Fx(z)-F 2 (z), zeTZ 3 , Rez<0. 



(4.5) 



Then these functions have an analytic continuation to a meromorphic function (denoted by £(z)) 
on the Riemann surface whose only pole is at the point at infinity on the first sheet. 

Proof. Let us first check that is analytic on the sheet TZk, k — 2,3. This reduces to showing 
the equality 

F %+ {z) + F 2 ,-(z) - Fx(z) = 0, z G flft\ [-ic,ic], 

which is a direct consequence of the variational condition (|2.15p . see also (|3. 31) . 

Next, we must show that the ^-functions are each others analytic continuation when passing 
through one of the cuts of the Riemann surface. For example, the fact that = £2,- along 
the cut R + n U*Li[ a i>&j'] reduces to showing the equality 



N 



F h+ (x) + F h -(x) - F 2 (x) - V'(x) + a = 0, 



x G 



n L)(°i. 6 i). 
3=1 



which is a direct consequence of the variational condition (|2.13[) . The other equalities are checked 
similarly, see also (|3.4p . □ 



It follows from Proposition 14.11 that the function £(z) is an algebraic function satisfying an 
equation of the third degree in £, known as the spectral curve: 



f+P2(z)e+Pi(z)^+Po(z) =0 
where po, p%, p 2 are polynomials. Here 

P2(z) = -Zi(z) - &(z) - &{z) = -V'(z) 



(4.6) 



(4.7) 
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is known, but the determination of the polynomials 

p (z) = 

cannot be done in general. We can only certify that (we use that V(z) is a polynomial of degree 
2d and FiO) = l/z + 0(l/z 3 ), F 2 (» = 1/(22) + 0(1/2) as 2 00) 

Pl (z) = V'{z)F 1 (z) - a 2 + O (z- 2 ) 

= YM. - a i + o (z 2d - 4 ) 

z ' 
p (z) = a 2 V'(z)~^l+O(z 2d - 3 ) 

as 2 — > 00. 

Example 4.2. For d = 1 we are in the quadratic case. If 




then pi(z) = 1 — a 2 , po(z) = a 2 z so that the spectral curve is 

f-z£, 2 + (l-a 2 )£ + a 2 z = 0. (4.8) 
This is known as Pastur's equation [31]. It plays an important role in [21 [HE]. 
Example 4.3. For d = 2we are in the quartic case. Let's take 

V(z) = y-\z 2 . 

Then pi(z) = z 2 + 0(1) and P o{z) — a 2 z 3 + O(z) as 2 — > 00. By symmetry we have that p\ is 
an even polynomial and po is an odd polynomial, so that pi and po are determined up to two 
constants. Thus for some a, G M, we have 

Pi(2) = 2 2 + a, Po(-z) = a2 z 3 + /32, 

and the spectral curve is 

£ 3 - (z 3 - tz) £ 2 + (z 2 + a)£ + a 2 2 3 + /?2 = 0. (4.9) 

This is McLaughlin's equation, named after K.T-R McLaughlin [3U] who derived it first for the 
case t = 0, see also [3J. We will analyze this case in more detail in Section [7] below. 

5 Proof of Theorem 12.41 

Recall the RH problem for Y(z) in ([2~2H|) - (|2T2T1) . In Subsections we will perform a 

Deift-Zhou steepest descent analysis of this RH problem. This will then lead to the proof of 
Theorem |2H in Section IET1 
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Figure 5: Lens around Supp(a — ^2) — *R in Case I. 



5.1 First transformation Y 1— > X 

In the first transformation we open up an unbounded lens around supp(cr — /12) which is bounded 
by a contour T = T + UT". We choose the lens so that it is symmetric under reflection with 
respect to both the real and the imaginary axis. The construction of the lens depends on whether 
we are in Case I or in one of the other cases (Case II or Case III). 

First assume that we are in Case I so that supp(<7 — /J2) = iBL Then we take T as in Figure[5] 
That is, we take q € (0,cln/2 + 1) and let r + be an unbounded Jordan curve in the right half- 
plane, symmetric in the real axis and crossing the real axis in q. Also T + is asymptotic to the 
half-ray arg z = ±8 for some 9 € (0, n/2). Then T~ is the reflection of T + in the imaginary axis. 
We orient T as shown in Figure [3] 

Next, consider the Cases II and III. Then we take T — T + U T~ as in Figure El The part 
of r + in the upper half plane is a Jordan curve going from 00 at an angle 9 € (0,7r/2) to the 
point ic. The other part of T + is its reflection in the real axis, and T~ is obtained from T + by 
reflection in the imaginary axis. We orient T as shown in Figure [5] 

The precise way to choose the contour T will be described further on. 

The contour T divides the complex plane into an inner and an outer part. By definition, 
we say that supp(er — /12) is inside the lens, while supp(/ii) is outside the lens. Note that our 
definitions are such that the outside of the lens is always on the left when traversing T according 
to its orientation. 

We define a new 3x3 matrix valued function X by 



X{z) 




Y\ :) I n I -r J . for Rez > outside the lens, 

(5.1) 

for Re 2; < outside the lens, 
for z inside the lens. 



Recall that q is used to denote the intersection of r + with the real axis (in Case I). In Cases 
II and III we put q = 0. 

Then X satisfies the following RH problem. 

(1) X is analytic in C \ (1 U [-ic, ic] U T). 
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Figure 6: Lens around Supp (a — ^2) in Cases II and III. 



(2) X satisfies the jump properties 

(1 e -n(V(x)-ax) q\ 
1 J for x > q, 

l) 

(1 e -".(V(x)+ax)\ 
1 for x < -q, 

1 J 

g— n(V(x)— ax) g — n(V(iE) + a:r)^ 

1 ) for - g < a; < q (in Case I), 

1 

on the real line, 

l\ 

= e- 2naz I for z e (-ic,ic) (in Case II, III), 

\0 -e 2naz 1 

on the interval (— ic, ic) oriented upwards, and 

/l 

X+(z)=X_(z) 1 - e - 2naz ] forzer+, 
\0 1 

/1 o^ 

X+(z) = 1 1 fazer 

\0 -e 2naz 1; 

on the contour I\ 
(c) As z -> 00 we have 



^" 
X(z) = (I + 0(l/z)) ( z-"/ 2 

Z-"/ 2 y 
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5.2 Second transformation X i— > T 

In the second transformation In>Twe use the minimizers /ii and /i 2 of the vector equilibrium 
problem, and the associated ^-functions 

9k{z) = J log(z - s)dfi k (s). (5.9) 

For k — 1 we choose the branch of the logarithm in (|5.9I) in the standard way, with a branch 
cut (-co, s] along the negative real axis. Thus g\{z) is defined and analytic in C\ (-co, fr/v] and 

gi(z) = log z + 0(1/ z), as z — > oo. 

For k — 2 the logarithm in (|5.9[) is taken with branch cut (—zoo, s] along the imaginary axis. 
Thus log(z — s) = log |z — s| + z arg(z — s) with — 7r/2 < arg(z — s) < 3it/2. Then 32(2) is defined 
and analytic in C \ iR, and 

.92(2) = ^ log 2 + 0(1), as z ^ 00, 

with the above branch of the logarithm. 

The behavior of the real and imaginary parts of the g-functions is described in the following 
lemma. 

Lemma 5.1. For x £ R \ {0} we have 

g lt+ (x)+ gi< -(x)=-2U^(x) 
gi,+(x) - gi-(x) = 2ni(i 1 ([x, +00)) 

0, x > 



g 2 (x) = -U^(x) + 



(5.10) 



7TZ/2, X < 0. 



For z€ jR\ {0} we /icwe 

S2,+(*) + ftJ,-(*) = -2C/^ 2 (z) +7ri/2 

27rz/i 2 ([z, ioo)), z e 



ff2,+ (z) - 02,- (*0 : 

ffi(z) = -^ I («) + 



ni/2 — 2az, z£[—ic,ic] (5-11) 

ni/2, z e zR+ 

-ni/2, z€zR_. 



Proof. The equalities follow immediately from the definitions, where we have to be careful with 
the choice of branches of the logarithm as discussed above. 

Only the equality in (|5.1ip for z £ [—ic, ic] needs some extra comment. We have that fj,% = a 
on [—ic, ic] and so if z £ [0, ic], 

Ha([0,z])=tr([0,z]) = ~\z\ = --iz, 

and therefore 

/i 2 ([z, zoo)) = M2QO, zoo)) - /i 2 ([0, z]) = - + -iz, (5-12) 

since \ii is symmetric in the origin and /i 2 (zR) = 1/2. The equality (I5.12[) is proved for z £ [—ic, 0] 
in a similar way, and thus 

92,+ (z) — g2,-(z) = 2ni^ 2 ([z, zoo)) = 7rz/2 — 2az, z £ [—ic, ic]. 

□ 
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Now we introduce the A-functions, which are defined as the following anti-derivatives of the 
£- functions (I4.5[) . Recall that I is the constant in the variational conditions (|2.13[) - (I2.14[) . 



\ 1 (z) = V(z) + l-g 1 (z), 
X 2 (z) = az + gi(z) - g 2 (z), 
\ 2 (z) = az + g 2 (z), 
A 3 (z) = -az + ni/2 + g 2 (z), 



z e n u 

z e TZ 2 , Rez > 0, 
zeTZ 2 , Rcz<0, 

z G 72-3, Rez > 0, 



(5.13) 



X 3 (z) — —az + ni/2 + gi(z) — g 2 (z), z G IZ3, Rez<0. 

Note that X 2 (z) and A 3 (z) are defined with a cut along the entire imaginary axis. In fact, from 
(j5TTTj) ([5TT3]) it follows that 



A3, + {z) - X 3 -(z) = 



0. 

— Til, 



z G [ic, ioo) 
z € (—ioo, - 



(5.14) 



and a similar formula holds for A2. 

We can now reformulate Lemma 15. II in terms of the A-functions. This leads to the following 
two lemmas. 



Lemma 5.2. We ho 



Re(A 3 ,+ (z) - 



- >>i,t( x ) 



= 0, x G M_|_ H supp(/ii), 

< 0, x G M + \ supp(/ii), 

= 0, x € M_ H supp(^i), 

< 0, x e M_ \ supp(/ii), 

= 0, z G iR \ (— ic, ic), 

< 0, z G (—ic, ic). 



(5.15) 



Proof. These are reformulations of the variational conditions (|2.13[) - (|2.16[) associated with the 
equilibrium problem, taking into account (|5.10p - (|5.13l) . The fact that we have strict inequalities 
follows from our assumption that the measure [i\ is regular. □ 

Lemma 5.3. We. have 



Al,4 


.(x) 


-Ai,- 


.(x) 


= -27ri^i([a:, +00)), 


X G 


^2,4 


\x) 


-A 2 ,_ 


.(x) 


= 2nifi 1 ([x, +00)), 


X Gl 


A3, -i 


.(x) 


-A 3) - 


_(x) 


= 2nifi 1 ([x, +00)), 


X G I 


\- 


v(z) 


-A 3) - 


_(*) 


= A 2 ,-(z) - A 3j+ (z) = 





(5.16) 



0. 



z G (—ic, ic), 
z G iK. 



2irifj, 2 ([z, ioo)) — ni/2 + 2az, 

The last expression is purely imaginary and its imaginary part is strictly increasing in terms of 
Im(z) as z G (—ioo, — ic) U (ic, ioo). 

Proof. This is a straightforward calculation using (|5.10l) - (|5.13l) . 



□ 



We define the new matrix valued function T as 



X(z) 



T(z) = diag (e-^,l,e-"^™ 

x diag ^ 

Then T satisfies the following RH problem 



e n(X 1 (z)-V(z)) e n(X 2 (z)-az) e n(X 3 (z)+az) 



(5.17) 
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(1) T is analytic onC\(RU [-ic, ic] U T). 

(2) The jumps for T are 



r e n(A 1 ,+ -Ai,_) e n(A 2 , + -Ai,_) q\ 

T+(x) = T-(x) | e »(A 2 , + -A 2 ,_) | fon>g ; (5.18) 

lj 

i e "(Ai,+ -Ai,_) q e n(A 3 , + -Ai,_)\ 

T+(x) = T_(x) ( 1 ) forx<-q, (5.19) 

e"( A3 .+- A3 '-), 

'\ e n(A 2 — Ai) e "(A 3 -Ai)\ 

T+(x) = T_(x) ( 1 j for - g < a; < 9 , (5.20) 

.0 1 



on the real axis, 



'10 

T+{z) = T-(z) \ 1 j for z€ (-ic,ic), (5.21) 

-1 e "(A 3 ,+ -A3,-) 



on the imaginary interval, and 



'10 

T+{z) = T-(z) ( 1 -e"( A3 - A2 ) ) forzer+, (5.22) 
v 1 

'i o o^ 

T+(z) = T_(z) I 1 0) forzer, (5.23) 



on the curve T. 
(3) As 2 — ► oo we have that 



_ e »(A2-A 3 ) i 



T{z) = I + 0{l/z). 



Note that for the (2, 3) and (3, 2) entries of the jump matrix in (|5.21[) . we used the last formula 
in (|5.16p . Also observe that T(z) does not have jumps on iM. \ [— ic, ic]. Indeed, in principle one 
should consider the jump relation 

(1 \ 

T + (z)=T-(z) e »(A2,+-A 2 ,_) o ze (-ioo,-ic)U(ic,i°°)- (5.24) 

\0 e n(A 3 ,+-A 3 ,-)y 

But from the jump relations in (|5 . 14[) we see that 

e „(A 3 ,+ (*)-A 3 ,_(z)) _ e o = lj z g ^ zoo ^ 

and also 

e r 1 (A 3 ,+ (z)-A 3 ,_(z)) = e -r 17 ri = ^ ^ g _ ^ 

by the fact that n is even. In a similar way one shows that e 71 ^ 2 ^^'^ 2 '-^' = 1, and so the 
jump matrix in (|5.24[) is just the identity matrix. 
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By virtue of (I5.15[) - (|5.16[) . we may rewrite the jumps on R as follows: 

x E supp(/ii), x > q, 



T+ = T_ 




1 ON 

i(Ai, + -Ai,_) q 







1, 



(5.25) 



and 



T+ = T_ 



s n(Ai I+ -Ai,_) q 



T+ = T_ 




-2-Kinotj „n 



(A 2l+ -Ai,_) q\ 



x e supp(//i). x < —q, 



, x € (bj,a j+1 ), x > q, 



(5.26) 



(5.27) 



T+ = T_ 



-2nina 



1 







where bo = — oo, aw+i = +oo and 

Oj = /ii([Oj+l, +Oo)), 



x e (bj,a j+ i), x < -q, 



j = 0,...,N. 



(5.28) 



(5.29) 



We would like the jump matrices on T to be exponentially close to the identity matrix as 
rt — > oo. From the above RH problem, we see that this is achieved provided T lies in the region 
where Re(A3 — A2) < if Re z > and Re(A3 — A2) > if Re z < 0. The fact that T can indeed 
be chosen in this way, follows by applying the Cauchy-Riemann equations to the last equality in 
Lemma 15.31 and using the last line in the statement of that lemma. 



5.3 Lenses around the intervals [aj,bj] 

Around each of the intervals [cij, bj] we open up a small lens to transform the oscillatory entries 
of the jump matrix into exponentially decaying entries. Since the non-trivial part of the jump 
matrix is locally of size 2x2 only, this can be done in the standard way [2] . 

More precisely, we take Jordan curves and TJ surrounding the interval [aj,bj] as in 
Figure The region between these curves is called the lens, and T'j' and TJ are the upper and 
lower lip of the lens, respectively. We choose them sufficiently close to the real axis so that 

Re(Ai(z) — A2(z)) < 0, for z inside the lens, z g [dj,bj], Re(z) > 0, (5.30) 

Re(Ai(z) — A3 (z)) < 0, for z inside the lens, z [o>j,bj], Re(z) < 0. (5.31) 

The fact that this is possible, follows from applying the Cauchy-Riemann equations to the first 
equation of ([STo]) , cf. [H] , 

The curves are chosen disjoint from each other, disjoint from the lens T in Section T5.ll and 
also disjoint from the imaginary axis. The only case where an intersection with the imaginary 
axis occurs is in Case III with j = (N + l)/2, and we assume that in this case, intersects iM. 
at a point ±izo strictly inside {—ic,ic). 
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Wc define 



1 O x 

TO) | -e™( Al ^- A2 W) 1 

1; 

/ 1 ON 

T(z) e "(^W-^(*)) 1 

V 1, 
SO) = <( / 1 0^ 

T(z) 10 

I _ e n(Ai(z)-A 3 (z)) ly 

/ 1 ON 

T(z) 10 

I e n(Ai(z)-A 3 (z)) Q i 

I ^0), 

Then 5 satisfies the following RH problem. 

(1) 5 is analytic in C \ (r U iR U T U \Jj Tj) 

(2) For a; e supp(^i) we have that 



1 ON 
= S-(x) (-100 

1, 

1 N 

S+(x) = S-(x) | 1 

,-1 0/ 

On the lips Tj of the lenses we have 



z in upper part lens, Re z > 0, 

z in lower part lens, Re z > 0, 

z in upper part lens, Re z < 0, 

z in lower part lens, Re z < 0, 
elsewhere. 



x e supp(/xi) n 



x e supp(/xi) n R_ 



1 ON 

S+O) - S_0) ( e"( Al - A2 ) 1 
1; 



(5.32) 



for^eTJ, Rez>0, 
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1 0^ 

S+(z) = S-(z)\ 10), for z erf, Rcz<0. 



e n(Ai-A 3 ) Q i 



In Case III we also have 

S+(z) = S-(z) for z € (— izQ,izo). (5.33) 

The jumps of S(z) on the other contours are the same as those for T(z). 
(3) As z — > oo, we have that 

S(z)=I + 0(l/z). 

From (|5.15[) and (I5.30[) - (|5.31[) , it can be checked that all the non-constant entries in the jump 
matrices for S(z) tend to as n — > oo, uniformly for z bounded away from the branch points dj, 
bj, ±ic. The only case that requires more explanation is the (3, 1) entry in the jump matrix in 
(|5.33[) . In that case, one can factorize 

e n(Ai(z)-A 3 ,_0)) = e n(Ai(*)-A 3 , + W) e n(A 3 , + (z)-A3,-W) ) 

and observe from (|5. 15|) and (15.31[) that for z £ [—izo,izo], the leftmost factor is uniformly 
bounded by 1 while the rightmost factor is uniformly exponentially decaying as n — > oo. 

5.4 Global parametrix 

The global parametrix we look for is a 3 x 3 matrix valued function M with jumps (obtained 
from the jumps of S by ignoring all entries which are exponentially small for n — > oo) 

(1) M is analytic in C \ ([ai, 6jv] U [— ic, ic\), 

(2) The jumps for M are 



x e supp(^i) n R+, (5.34) 



x € supp(/xi) nl_, (5.35) 

U-(.c) = M-(.r) I e 2™a 3 oj x € (p j:CLj+1 ) n iR +! (5.36) 
1/ 

2Triria.j q q N 

M + (x)=M-(x)\ 1 ), x G (bj, aj+i) n R_, (5.37) 

e 27rmQ ^ 

on the real line, and 






1 
































1 














'\ 0\ 

M+(z) = M_(z) | 1 , z e (-ic,ic). (5.38) 
,0 -1 0/ 
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Figure 8: Canonical homology basis in Case I (with N = 4). 




Figure 9: Canonical homology basis in Case II (with N = 4). 



(3) At infinity we have 

M(z) = I + 0(1/ z), as z -t oo, 

(4) M has at most fourth-root singularities at the branch points oi, b\, . . . , Qn, 6 at, ±ic. 

We can solve this problem with the help of meromorphic differentials on the Riemann surface. 
Such a construction was first used in [3T] and later developed further in [55] . 

To the Riemann surface we associate a canonical homology basis {A±, . . . , A g , B±, . . . , B g } 
where g is the genus. The details of the construction depend on whether we are in Case I, II or 
III, see Figures EHini 

For brevity, we give a detailed description only for Cases II and III. Then the genus is g = 
N — 1. The cycles Bj are on the first sheet and Bj encircles [oi, bj] once in the counterclockwise 
direction. The cycles Aj are partly in the upper half-plane on the first sheet and partly in the 
lower half-plane on the second or third sheet. Aj passes through [<Xj,6j] and [dj+i, 

The anti-holomorphic involution (j) is defined by mapping z to z on the same sheet. The 
fixed point set of (f> is the disjoint union of N = g + 1 closed curves (Jf=o on ^ ne Rie maim 
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Figure 10: Canonical homology basis in Case III (with N = 3). 



surface. Here Sj is nomotopic to Aj as a closed curve, j = l,...,g, while So is the unbounded 
component. 

Lemma 5.4. If Pj G £j for j = 1, . . ,,g, then the divisor 



is non-special. 

Proof. A detailed proof of this theorem in a similar setting will be given in [2T] . Here we only 
outline some of the key steps. We use u : 1Z —> C 9 /L to denote the Abel map, mapping the 
Riemann surface 1Z to <C 9 /L, where L is the lattice defined from the columns of the Riemann 
period matrix. We also use 9 : C 9 — > C 9 to denote the corresponding Riemann theta function on 
C s , and K denotes the vector of Riemann constants, see e.g. [25] . 

The proof is based on the following result which can be found e.g. in [T5J Theorem 2.4.2]. 
The divisor D is non-special if and only if 9(u(P) — u(D) — K) does not vanish identically for P 
on the Riemann surface. 

The next step is to show that for s € M 9 we have 9{s) > 0. This relies on the fact that the 
Riemann surface 1Z is the Schottky double of a bordered Riemann surface, see [Ml Corollary 
6.13] and [31]. 

Using the antiholomorphic involution <j) it can be shown that the Riemann period matrix is 
purely imaginary. Finally, one then shows that u(D) + K has a real representative modulo the 
lattice L [21] . By taking into account the results in the last two paragraphs, the desired result 
then follows. □ 

We now basically follow [25]. Given (Pi,...,P g ) with Pj £ Sj, we define a meromorphic 
differential up so that 

• cop has simple poles in a\, b\, . . . , ajv, &jv, ±«c, Px, . . . , P g , 002 and 003 with residues 




i=i 



Res(o;p,aj) — Kes(u>p,bj) = — |, 
Res(wp,±ic) = — 7j, 



j = l,...,N, 



Bss(u} P ,Pj) = 1, j = l,...,g 
Res(wp, 002) = Res(ajp, 003) = 1, 
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• uip is holomorphic elsewhere, 



• lup has vanishing A-periods: 



L 



LUp = 0, 



3 = 1, 



The total sum of the residues is 0. These conditions determine uip uniquely. The B-periods 
are purely imaginary, they are well defined modulo 2-7riZ, and the mapping 



Let ujp be the corresponding meromorphic differential. 

We take the base point Pq = ooi and define three functions v\(z), V2(z) and v^(z) of a 
complex variable z as follows. We have 



\JPo J 

where z is considered as a point on the /cth sheet of 1Z, and where the path of integration is as 
follows 

• for k = 1, the path of integration is on the first sheet and does not intersect the real line, 

• for k = 2, the path of integration is on the first and second sheets; for Imz > (Imz < 0), 
the path starts in the lower (upper) half-plane of the first sheet, crosses one of the cuts 
[aj,6j] in M + to the upper (lower) half-plane of the second sheet and stays there until it 
ends at z\ 

• for k = 3, the path of integration is on the first and third sheets; for Imz > (Imz < 0), 
the path starts in the lower (upper) half-plane of the first sheet, crosses one of the cuts 
[dj, bj] in K_ to the upper (lower) half-plane of the third sheet and stays there until it ends 
at z. 

Then the vector (vi, i>2, W3) is well-defined and analytic in C \ (M U iR) with jumps 






mod Z. 




{vi,V2,V 3 )+ = {v!,V2,V3)-J-> 



v 



where the jump matrices on the real line are 







2G 




for x < ai and x > b 



N, 



for bj < x < a, + i, x > 0, 



for bj < x < Oj+i, a; < 0, 



while on the imaginary axis they are 



'\ Cf> 

Ju(z) = I 1 I , for z G (— ioo, —ic) U (ic, ioo), 




J v (z) — \ 1 , for z e (— ic, ic). 



Also as z — » oo, we have 

Wl (z) = l + 0(l/z), wa(«) = (9(1/2), v 3 (z) = 0(l/z), 

and all Vj functions have fourth-root singular behavior at the branch points 01, b\, ... , aAr, 6jv, 
±ic. 

Now define the first row (M11, M12, M13) as follows 

(M u ,Mi2,Mi 3 ) = («i,u 2 ,W3) if Imz > 0, 

= (vi,— V2,— va) iflmz<0. 

This gives the correct jumps for M. The other rows of M can be constructed in a similar way, 
or by a simple transformation of the first row . 

5.5 Local parametrices 

In the regular case (the one we are considering) we construct local parametrices out of Airy 
functions near each of the branch points. We denote the local parametrices by P. The local 
parametrices match with the global parametrix on the boundary of a small circle around the 
branch points. Since the non-trivial part of the jump matrix is locally of size 2x2 only, and 
since in the regular case we have square root behavior near all the branch points (Lemma 12. 3[) , 
this construction can be done in the usual way [14]. We omit the details. 



5.6 Final RH problem 

We define the final RH matrix 

. _ J T(z)M(z)^ 1 away from the branch points, 

1 T(z)P(z)^ 1 near the branch points. 

Then R(z) has jump matrices tending to the identity matrix as n — > 00 and it is normalized in 
the sense that R(z) = I + 0(1/ z) as z — > 00. Standard arguments then show that R(z) itself 
also tends to the identity matrix, uniformly for z on compact subsets of C [14] . 
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5.7 Proof of Theorem EH 

Having performed the steepest descent analysis of the RH problem, we can now prove Theo- 
rem 12.41 The proof follows the same pattern as in the papers [2j |8l [9j [18] . 

First assume that (a,j,bj) with x,y > 0. We will transform (|2.22p under the series of 

transformations Y m- X n- T i->- S. By virtue of (15. ip we have 

K n (x,y) = —± -(0 w x {y) 0) X+^X+fjc) ( 

2tti(x-?/) ^ 

Using (|5.17l) this becomes 

c »(v(.)-vto)) / e -*i.+<« > 

K n (x,y) = — — (0 e nX *>+M 0) T + ( 2 /)- 1 T + (x) 

2m{x-y) \ 

From ([5321) we get 



e n(V(x)-V(y)) 



' e 



KJx,y) = — — — (- e " A ^+(f) e " A2 .+^ 0) S' + (?/)- 1 S' + (x) e -" A2 .+W . (5.39) 

2m(x — y) 1 1 

Now it follows by standard arguments (e.g. Section 9]) that 

S^ 1 (y)S(x) = I + 0(x - y), as y ->• a; 



uniformly in n. Inserting this into (15 -39[) and setting h(x) = V(x) — Re(Ai(x)) yields 

K n (x,y) = e nW.)-%» ^("MAi,+ (»)-A 1>+ (g))) + ^ ^ 

Tr(x-y) 

uniformly in n. By letting y x and using l'Hopital's rule we find 

7T 

as n — > oo. From (14.51) and the Stieltjes inversion principle we conclude that 

lim —K n (x,x) = df*i( x ) _ 
7i ^ oo n ax 

This proves (|2.23l) for x € (a,j, bj), x > 0. The proof for the other values of x £ K is similar, or 
can be obtained from symmetry considerations. 



6 Phase transitions: General discussion 

Recall the Cases I, II and III describing the structure of the equilibrium measures in Section I3~3"l 
Intuitively one expects the following possible behavior in terms of the parameter a. For large a 
we are in Case I. The measure fi\ is supported on two (or more) disjoint intervals with a gap 
around 0. The constraint a is not active. 

When a decreases the gap around shrinks. Then one of two things could happen. It could 
happen that for a certain value of a the constraint becomes active, while the gap in the support 
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of [i\ around is still there. Then we are in Case II. Then if a further decreases the gap may close 
or not. The latter depends on whether in the unitary matrix model with potential V (without 
external source) is in the support or not. If the support closes then we are in Case III. 

The other situation that could happen is that the constraint a remains inactive all the way 
until for a certain value of a the gap in the support of jiii is closed. Then if a further decreases 
the constraint becomes active. The transition is then from Case I to Case III without passing 
through the Case II. This is precisely what happens in the quadratic case V(x) — \x 2 . More 
generally, one expects this kind of behavior when the potential V(x) is convex, or 'nearly' convex. 

For those values of a for which a transition between one of the Cases I, II, III to another 
takes places, one expects that the local eigenvalue correlations near the origin are described 
by special functions related to ODE's. For typical cases one expects such special functions as 
Pearcey integrals [TJ |H and the Hastings-McLeod solution to the Painleve II equation [5J [T3] . 
However, our model allows for new kinds of critical and multi-critical behavior as well, but it 
remains an open problem to describe these new critical phenomena. 

In Section [7] we will illustrate the above considerations in detail for the case of a quartic 
potential. See in particular Figure [TT1 

7 A case study: The quartic potential 

Let us investigate the case of a quartic potential 

V(z) = -z 4 - -z 2 
W 4 2 

and the associated McLaughlin equation (|4.9I) : 

£ 3 - (z 3 - tz) £ 2 + (z 2 + a)£ + a 2 z 3 + f3z = 0. (7.1) 

The discriminant of the McLaughlin equation (w.r.t. £) is a polynomial D\2{z) of degree 12 in 
z. We calculated it with Maple, but it is too long and not too interesting to reproduce it here in 
full. The first terms are 

D 12 (z) = -4a 3 + (a 2 t 2 + mafit - 12a 2 - 27 /3 2 )z 2 + 0(z 4 ) as z -> 0. (7.2) 

The branch points of the Riemann surface are among the zeros of D\2{z). There are other 
zeros, and they should come with higher multiplicities. 

For general a and (3 the McLaughlin equation has genus 4 (according to Maple). The special 
choices for a and /3 that are relevant to us will lead to a reduction of the genus. The genus can 
be at most one, as the following lemma shows. 

A similar result occurs in [501 Prop. 5.2.5], but the proof given there is incorrect. Here we 
give a self-contained proof which may be used for the situation in [2U] as well. The proof uses 
an idea due to Lun Zhang (personal communication). 

Lemma 7.1. Assume that x i— > V(y/x) is convex for x > 0. Then the support of fj,± is either 
one interval (in case € supp(/ii)J or a disjoint union of two intervals (in case ^ supp(/zi)^ 
and the measure \i\ can have singular behavior only at zero. 

Proof. By fixing fi2 < o in the energy functional (|2.1[) . we see that [i\ is the unique minimizer 
for the energy functional 

/oo 
(V{x)-a\x\-U f " 2 (x))d f i(x) 
-OO 



29 



among probability measures [i on R. Since the external field is symmetric, it follows from [351 
Theorem IV. 1.10(f)] that d^i(t) = djii(t 2 )/2 where j2\ is the unique minimizer of the energy 
functional 

/>oo 

/(M) + 2 / {V{^)-a^-U^{^))d^{x) (7.3) 
Jo 

over all probability measures /i on [0, oo). 

We are going to show that the external field V(y / x) — a^Jx — U^ 2 (y/x) in (|7.3p is convex for 
x > 0. Since 

U^{y/X) = ~ J \og{x + \z\ 2 )d m {z), 

we have 

which due to the constraint /j,2 < cr can be bounded by 

d 2 , ,, , /— ,\ a /" 1 , , , a 

^ (u^) < Tv I 1*1 = ^ 

= ^ (-«v^) • 

It follows that —a^/x — U^ 2 (\/x) is convex. Due to the assumption in the lemma on the convexity 
of V(</x), it then follows that (V(^/x) — a^Jx — C/ AI2 ( v / x)) is indeed convex for x > 0. 

From the convexity it follows that jli is supported on one interval. Then /ii is supported 
on either one or two intervals, depending on whether belongs to the support of jli or not. 
The convexity also implies that the density of /ti does not vanish in the interior of its support, 
and has square root behavior at its non-zero endpoint(s), see e.g. |12[ Lemma 3.5]. The same 
properties then apply to fX\, □ 

The three cases in Section [3~31 now specialize as follows. 

Case I: In the first case there are four real branch points ±&i, ±62 (with b\ > 62 > 0) and no 
other branch points. The remaining zeros of the discriminant (|7.2H come as four double zeros. 
The Riemann surface has genus zero. 

Case II: In the second case there are four real branch points ±61, ±62 (with 61 > 62 > 0) and 
two purely imaginary branch points ±zc (with c > 0). The remaining zeros of the discriminant 
come in the form of a six- fold zero at z = 0, see Lemma 17.21 below. The Riemann surface has 
genus one. 

Case III: In the third case there are two real branch points at ±61 (with bi > 0) and two purely 
imaginary branch points at ±zc (with c > 0). There are four double zeros of the discriminant. 
The Riemann surface has genus zero. 

It remains to investigate in more detail the Cases I, II and III. More precisely, for any i e 1 
and a > we want to determine which of the three cases applies. We also want to find the curves 
in the (t, a) plane where a transition from one case to another takes place. 
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The genus one region 

Let us first investigate Case II. This is the case of genus 1. It turns out that in this case, the 
parameters a and f3 in the McLaughlin equation (|7.1|) take on a particularly simple form: they 
are both equal to zero. This is the content of the next lemma. 

Lemma 7.2. (Case II.) Let f 6 1 and a > be such that the McLaughlin equation (|7.1|) is of 
genus 1. Then one has that a = j3 = 0. Moreover, the discriminant of the McLaughlin equation 
can then be factored as 

D 12 (z) = z 6 D 6 (z) (7.4) 

where Dq is the degree six polynomial 

D 6 (z) = Aa 2 z 6 + (1 - Vita 1 ) z 4 + ((12< 2 - 18)a 2 - 2f) z 2 

+ (-27a 4 + (18*-4< 3 )a 2 -A + t 2 ) . (7.5) 

Proof. From the general descriptions above we know that in the genus 1 case, the McLaughlin 
equation has six simple branch points ±6i, ±62, and ±ic (61,62,0 > 0), which are simple zeros 
of the discriminant. The remaining zeros of the discriminant should come as three double zeros, 
possibly coalescing. By symmetry is a double zero, and this forces a = 0, cf. (|7.2[) . 

If we substitute a = in (|7.1|) and calculate the discriminant with respect to £ we obtain the 
12th degree polynomial 

4aV 2 - (-4/3 + 12a 2 * - l)z 10 - (18a 2 - 12< 2 a 2 + 2t + I2tfi)z 8 
- (4 - 12t 2 f3 - t 2 + 18/? - 18a 2 t + At 3 a 2 + 27a 4 )z 6 - (54a 2 /3 - 18*/3 + 4* 3 ,3)z 4 - 27/3 2 z 2 . (7.6) 

We can take out the factor z 2 in (|7.6[) . corresponding to the double zero at z = 0. The remaining 
factor is of degree 10 and we consider it as a 5th degree polynomial D§{y) my — z 2 . We already 
know that this polynomial has three simple roots 6 2 , b\, — c 2 and one double root, call it d. 
Since by (|7.6[1 the sum of the roots equals (—4/3 + 12a 2 * — l)/(4a 2 ) € R, it follows that d is real. 
Moreover, from (I7.6[) we see that the product of the roots equals 27/3 2 /(4a 2 ) > 0. Since this 
product can also be written as — 6 2 6|c 2 d 2 < 0, it then follows that /3 = d = 0. 

In conclusion, we have shown now that a = (3 = 0. Inserting this in the McLaughlin equation 
and computing its discriminant by a direct calculation then leads to (|7.4l) - (|7.5p . □ 

Note that the factor z e in (|7.4I) corresponds to the six- fold zero at z = Q, while the zeros of 
Dq(z) should yield the branch points ±61, ±62 and ±zc. In particular, four of these zeros should 
be real and the other two purely imaginary. The next lemma describes when this happens. 

Lemma 7.3. The polynomial (|7.5[) has four real and two purely imaginary zeros precisely for 
those (t, a) € Kx R + lying in the open region T> bounded by the points (ti, ax) = (2, 0), fa, a 2 ) — 
l/v27), (£3,(13) = (+oo,0), the straight line segment between (ti,Oi) and (£3,03), and the 
two branches of the curve 

As(0) = -27a 4 + (18* - 4i 3 )a 2 -A + t 2 = (7.7) 

that connect the point (£2,02) with (t±,ai) and (£3,03), respectively. The region T> is shown in 
the bottom right part of Figure llll 

Proof. Rewrite Dq(z) as a cubic polynomial in the variable y — z 2 : 
4a 2 y 3 + (l - 12ia 2 ) y 2 + ((12t 2 - 18)a 2 - 2t) y 

+ (-27a 4 + (18* - 4* 3 )a 2 -A + t 2 ) . 
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Figure 11: Phase diagram for the quartic potential jx 4 — |x 2 . The dark curve represents the 
Painleve II transition. The dotted curve represents the Pearcey transition. 

We are looking for the values of (t, »)elx R + for which this cubic polynomial has two strictly 
positive and one strictly negative zero. This is a standard routine whose description we omit. □ 

The two previous lemmas show that the genus of the McLaughlin equation can only be 1 if 
(t, a) lies in the region T>. Outside T> the genus must necessarily be zero. 

It remains to show that inside T> the genus is exactly 1 (and not 0). This is taken care of by 
the next lemma. 

Lemma 7.4. (a) Inside T> the genus is either identically 1 or identically 0. 

(b) There exists at least one point (t, a) in T> for which the genus is 1. 

Proof. For Part (a) we invoke the following continuity result which follows from general potential 
theory: The parameters a and (3 in the McLaughlin equation depend continuously oniGt and 
a > 0. Using this, one sees that the discriminant of the McLaughlin equation also depends 
continuously on t and a and therefore also the branch points since these are simple roots of the 
discriminant. 

Now let D\ C T> be the region formed by those (t,a) £ T> for which the genus is 1. We 
show that T>i is both open and closed in V. To show that it is open, let (t, a) £ T>\. Then the 
discriminant has six simple zeros and by continuity the same holds in an open neighborhood of 
(t, a). To show that T>\ is closed in T>, we take a sequence of points (ife, flfe) € Pi, = 1,2,..., 
which converge to a limit point (t, a) £ V. By Lemma P7T21 we have a = ft = for each (tfe, a&) so 
by continuity the same must hold for the limit point it, a). But then Lemma 17.31 shows that the 
discriminant has six distinct simple zeros, which implies that the genus is 1. Hence (t,a) £ T>\. 

For Part (b), we only outline a proof. The idea is to show that for any fixed t > 2, we have 
(t, a) £ T>\ for all a small enough. This relies on the fact that for t > 2 the eigenvalues in the 
unitary matrix model with potential jx 4 — ^x 2 (without external source) are supported on two 
intervals [UE]. The claim then follows from a continuity argument for a — > 0; we do not go into 
the details. 

An alternative approach to prove Part (b) would be to pick a numerical point (i, a) £ T> and 
show by direct means (using the McLaughlin equation with a = (3 = 0) that this algebraic curve 
makes the RH steepest descent analysis work, in a similar vein as in [51 [3] ■ □ 
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We summarize our findings with the following 



Proposition 7.5. Let T> be the region defined in Lemma 17.31 For (t,a) 6 T> the McLaughlin 
equation has genus one and we have a = /3 = 0. For (t,a) £ (R x R + ) \ 2? t/ie genus is zero. 

The genus zero region 

Now we focus on the genus zero region (t, a) E (R x R + ) \ 27. A useful representation for 
the parameters a and (3 in this region is given in [3]. It is shown there that the following 
parametrization holds 

(l-u)(c 4 + 2c 4 u-u 2 ) 2 
a = — 7.8 

p = J3c 4 -u)( c 4 + 2 c 4 u-u 2 ) 2 

C 4 M 2 

where 

u:=2c 4 + ac (7.10) 

and where c is a root of the equation 

2c 6 - 2<c 4 + 3c 2 - ac 3 - tac - a 2 = 0. (7.11) 

According to [3J, c is actually the largest positive root of equation (I7.11|) . but we will not need 
this in what follows. 

Let us seek the values of (t, a) £ (R x R + ) \ V for which a phase transition occurs between 
the Cases I, II, III. For such (t, a) the discriminant of the McLaughlin equation should have a 
zero at z = 0. By virtue of (|7.2p this implies a = 0. From (|7.8[) this implies that either u = 1, 
or c 4 + 2c 4 m — u 2 = 0. Let us first consider the case u = 1. Then we have from (|7.10[) that 

2c 4 + ac=l. (7.12) 

Thus the equations (|7.1ip - (|7.12l) have a common root c. In other words, the resultant of these 
two equations with respect to the variable c should be zero. Computing this resultant with Maple 
yields the following condition on (t,a): 

54a 4 + (72t - t 3 )a 2 - (t 4 - I6t 2 + 64) = 0. (7.13) 

Next we consider the case c 4 + 2c 4 m — u 2 = 0. From (|7.8p ~ (|7.9p this implies that a — f3 = 
and we know from earlier considerations (or from a similar resultant calculation as above) that 
this is only possible if (t, a) is such that (|7.7p holds. The phase transition on this curve will be 
discussed in the next section. The phase transition on the curve (|7.13l) will be discussed in the 
section thereafter. 

Painleve II transition 

At the two curved boundaries of the region T> in Figure II H we have a transition from genus to 
genus 1. Recall that these boundaries are described by the relevant branches of the equation 

D 6 {0) = -27a 4 + (18* - 4< 3 )a 2 - 4 + t 2 = 0. 
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More precisely, these branches are given by 



Ai = a 2 



An 



t 2 

3 ~ 27 
t 2 

3 ~ 27 



t 3 _ ( t 2 _ 3)3/2 
t 3 + (t 2 _ 3) 3/2 



t > \/3 
V3 < t < 2. 



We have A\ > for every i > whereas A 2 > only for \/3 < t < 2 and < A 2 < Ax for 
these i- values. See Figure [TTJ 

On the above curves we have a transition from genus to genus 1, and we expect that the 
phase transition is of Painleve II type [5J [T21 HB] . More precisely, we expect that the following 
happens. If one lets a decrease towards the curve a 2 = A\ (with t > \/3), then the constraint a 
on the imaginary axis becomes active, and we have a transition from Case I to Case II. If one 
further decreases a towards the curve a 2 = A 2 (V3 < t < 2), then the gap in the support of fix 
closes and hence we have a transition from Case II to Case III. 

On the curve a 2 = A 2 the phase transition involves the eigenvalue measure \i\ and therefore we 
expect Painleve II behavior in the local eigenvalue correlations at the origin [51 H2] . On the curve 
a 2 = Ax, however, the phase transition takes place on the 'non-physical' sheets of the Riemann 
surface and therefore it is not felt in the eigenvalue statistics. But then we expect Painleve II 
behavior in the recurrence coefficients for the associated multiple orthogonal polynomials, as 
in PI]. 

Note that for t > 2 the transition at a 2 — A 2 does not occur. This is consistent with the fact 



that for t > 2 the eigenvalues in the unitary matrix model with potential . 
external source) are supported on two intervals [H[S], as mentioned before. 

Pearcey transition 

Now we consider the equation (|7.13l) . 

54a 4 + (72t - t 3 )a 2 - (i 4 - 16t 2 + 64) = 0. 
This equation has the two solutions 

A 3 =a 2 = 

A 4 =a 2 = 



ix 4 



(without 



2t 


1 - 


3 H 


" 108 . 


2t 


1 - 


3^ 


" 108 . 



t 3 + (t 2 + 24) 3/2 
t 3 - (t 2 + 24) 3/2 



The branch A4 is negative, and so is irrelevant for us. The other branch A3 is positive and 
we expect that for t < V3 a phase transition of the Pearcey type [TJ |H1 1331 121] takes place for 
a 2 = A 3 . See Figure [TT1 
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